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In this paper we present the new concept of description of random processes in complex systems 
with the discrete time. It involves the description of kinetics of discrete processes by means of 
the chain of finite-difference non-Markov equations for time correlation functions (TCF). We have 
introduced the dynamic (time dependent) information Shannon entropy Si(t) where i=0,l,2,3,... 
as an information measure of stochastic dynamics of time correlation (i — 0) and time memory 
(i=l,2,3,...). The set of functions Si(t) constitute the quantitative measure of time correlation 
disorder (i — 0) and time memory disorder (i=l,2,3,...) in complex system. The theory developed 
started from the careful analysis of time correlation involving dynamics of vectors set of various 
chaotic states. We examine two stochastic processes involving the creation and annihilation of 
time correlation (or time memory) in details. We carry out the analysis of vectors' dynamics 
employing finite-difference equations for random variables and the evolution operator describing 
their natural motion. The existence of TCF results in the construction of the set of projection 
operators by the usage of scalar product operation. Harnessing the infinite set of orthogonal dynamic 
random variables on a basis of Gram-Shmidt orthogonalization procedure tends to creation of infinite 
chain of finite-difference non-Markov kinetic equations for discrete TCF and memory functions 
(MF). The solution of the equations above thereof brings to the recurrence relations between the 
TCF and MF of senior and junior orders. This offers the new opportunities for detecting the 
frequency spectra of power of entropy function Si(t) for time correlation (i = 0) and time memory 
(i=l,2,3,...). The results obtained offer considerable scope for attack on stochastic dynamics of 
discrete random processes in a complex systems. Application of this technique on the analysis 
of stochastic dynamics of RR-intervals from human ECG's shows convincing evidence for a non- 
Markovian phenomemena associated with a peculiarities in short and long-range scaling. This 
method may be of use in distinguishing healthy from pathologic data sets based in differences in 
these non-Markovian properties. 
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I. INTRODUCTION 

Manifold methods are successfully used in statistical physics for the description of distinctive characteristics of 
chaotic dynamics of complex systems [|lj- |^6j . Nevertheless, three vexing features which are difficult to yield a detailed 
and strict analysis are available in complex systems. Among them: nonstationarity, nonlinearity, and nonequlibrium 
phenomena. Furthermore, the significant peculiarities of complex systems are directly related to the discretness of 
time of object / subject registration response ||, Q, |Q, |l7|] , Non-Markov and long-range statistical memory 
effects also play the leading part in the complex systems behaviour 0- j|, |l7t| , |p7j , fl33f . 

However, the discretness of time while considering the complex systems has not been taken into account until now, 
although it is discretness that is the most commonly encountered feature of real objects / subjects. On the other hand 
the memory and time long-ranging effects are paramount. As a rule the state developed is complicated by the fact 
that the real complex systems are of nonphysical nature. Therefore, the direct methods of statistical physics derived 
from Hamiltonian formalism, exact equations of motion and Liouville's equation are not applicable in this case to its 
theoretical analysis. Meanwhile the real existence of complex systems in time and space generates a reliable evristic 
basis for the modeling in terms of the time discretness, memory and time long-range effects. 

The present article is dedicated to statistical consideration of a discretization in temporary changes of complex 
systems of a substantial nature on the basis of the first principles. In Section 2 we briefly outline general definitions 
and proposals used to form the stochastic dynamics of discrete time sequences, and in Section 3 we suggest the 
geometrical presentation of stochastic dynamics of time correlation. Introduction of projection operators, splitting of 
equation of states vectors and matrix presentation of Liouville's quasioperator for the statistical description of random 



processes with discrete time are reported in Section 4, and introduction of the set of orthogonal random variables as 
well as construction of infinite chain of finite-difference non-Markov kinetic equations for discrete TCF are framed 
in Section 5. A pseudohydrodynamic description of random processes is provided in Section 6, where the relative 
merits of this approach are set forth. In Section 7 we define Shannon dynamical (time dependent) entropy for time 
correlation and time memory in complex systems. Aplication of technique on the analysis of stochastic dynamics of 
RR-intervals from human ECG's are discussed in Sec. 8. In Section 9 we present the discussion and conclusions of 
the results obtained and possible opportunities for the experimental data processing. 

II. BASIC ASSUMPTIONS AND DEFINITIONS 

Following Gaspard [15] we consider a random process such as a sequence of random variables defined at successive 
times. We shall denote the random variable by 

X = {x(T), x(T + t),x(T + 2t),---, x(T + fcr), ■ • • , x(T + tN - r)}, (1) 

which corresponds to signal during the time period t = (N — 1)t where r is time interval of signal discretization. The 
mean value < X >, fluctuations 6xj, absolute (a 2 ) and relative (S 2 ) dispersion for a set of random variables (|lj) can 
be easily found by 



<x>=^j2< T+ ^ ( 2 ) 
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Xj = x(T + jr), Sxj = Xj— < X >, (3) 
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Abovementioned values determine the statical (independent from time) properties of the system considered. The 
normalized time correlation function (TCF) [1-3] , [7-9] depending on discrete time t = tot, N — 1 > m > 1 can be 
conveniently used for the analysis of dynamic properties of complex systems 

j N—l—m 

"(*)= (N _ m)a 2 E Sx(T + jr)Sx(T + (j + mr)). (6) 

^ ' 3=0 

TCF usage means that developed method is just for complex systems, when correlation function exist. In forthcoming 
papers we intend to apply developed method for discrete random processes analysis in complex systems in practical 
psychology, cardiology (for the development of diagnosis method of cardiovascular diseases) , financial and ecological 
systems, seismic phenomena, etc. The properties of TCFs a(t) are easily determined by Eq. (^|) 

\\ma(t) = 1, lim ait) = 0. (7) 

We have to recognize that the second property in Eq. (]?]) is not always satistifed for the real systems even with 
arbitrary big values of time t or number [N — 1) = t/r. Taken into account fact that the process is discrete, we must 
rearrange all standard operation of differentiation and integration [34,35] 



The first derivative on the right is recorded in Eq.(||). The second derivative on the right is also derived easily 
<fx(t) Ax ( Ax" 



,11- • At V At ) = T ~ {[x{t + 2t) " X{t + T)] ~ [x{t + T) ~ X{t)]} = {x{t + 2r) " 2x{t + T) + X{t)} - (9) 

Now let us proceed to the description of the dynamics of the process. For real systems values Xj = x(T + jr) and 
Sxj — 5x(T + jr) result from the experimental data. Thus we can introduce in Shannon's manner [17] the evolution 
operator U(T + t%, T + 1\) in as follows (t2 > ti) 

x(T + fa) = U(T + t 2 ,T + n)x(T + t x ). (10) 



For brevity let us present Eq.(10) in the form 

x(j) = U(j,k)x(k), j > k, j,k = 0,1,2, - ■■ ,N — 1. (11) 

Time operator of one step shift r along a discrete trajectory is conviniently considered by means of two nearest values 
x(t + r) and x{t) 

x{t + r) =U(t + r,t)x(t). (12) 
Owing to Eqs.([To|)-(|T2[) a formal equation of motion is derivable for any x £ (xq, X\, x<i, ■ ■ ■ , xn-i) 



§ - ^ = r- l {x{t + r) - x(t)} = r-i{U(t + r, t) - l}x(t). 



(13) 



Let us consider Eqn.(fL3j) in terms of x 



A^f) = x 3+1 {t + T)-x 3 (t) = T _ 1{u{t + ^ t) _ 1}jBj(t) 



and then introduce a Liouville's quasioperator L as follows 

Ax(t) 



dt At 



iL(t,T)x{t),L{t,T) = {iT)- L [U(t + r,t) - 1]. (14) 



Now in accordance with Refs. (jf5[, |T(| let us present a set of values of random variables Sxj = Sx(T + jr),j — 
0,1, - ■ ■ , N — lasa fc-component vector of system state 

Ag(0) = {6x ,6x 1 ,6x 2 , • • • , Sxk-t) = (Sx(T), Sx(T + r), • • • , Sx(T + (k - l)r). (15) 

Now we can introduce the scalar product operation 

fc-i 

< A B >=^A j B j (16) 
j=0 

with or whithout indication of obvious time dependence of vectors A and B, respectively, in the set of vectors A°,(0) 
and A™ +fc (t) where t — mr and 

■^■m+fe W = {^ X nn 3 x m+l, " " j & x m+k-l} 

= {Sx{T + mr), Sx(T + (m + l)r), Sx(T + (m + 2)r), • • • , fe(T + (to + fc - l)r}. (17) 

A number k < N — 1 determines the vectors' dimension. The functions (0), (^|) can be expressed in terms of scalar 
product ( |I6| ) 

- 2 = ^<A° r .A^>=7V- 1 {AU 2 , 
r2 TV- 1 < A° • A° > 



A™ +fc (t) = U(T + mr, T)A°(0) = {U(T + mr, T + (m - l)r) 
E7(T + (m - l)r, T + (m - 2)r) • • ■ J7(T + r, T}A£(0). 



(18) 



The normalized TCF in Eq. (§) can be rewritten in a more compact form by means of Eqs.(|l^)-(|T8]) ( t = mr is 
discrete time here) 



a(t) 



<K-K l+k > <A2(o)-A™ +fc (t)> 



<A° fc -A°> 



< A°(0) 2 > 



Replacement of Eqn.(||) by Eqn.(|l9|) is true if the numbers k < N — 1 satisties the condition 

fc-i fc-i 
a 2 = k- 1 V Sx 2 or cr 2 = lim AT 1 V Sx 2 . 

The condition of quasi-stationarity of processes under consideration 



da(T,t) 



dT 



da(t) 



dt 



(19) 



(20) 



(21) 



serves the other criterion of validity for such replacement. TCF a(T,t) at Eqn.(|2lj) is viewed on a time scale (point 
T) at the distance t from the zero point . 

Such vectors' notion is very helpful for the analysis of dynamics of random processes by means of finite-difference 
kinetic equations of non-Markov type. 



III. GEOMETRICAL NOTION OF STOCHASTIC DYNAMICS OF TIME CORRELATION 

First of all let us consider the projection operation in the set of vectors for different system states. It is easy to 
introduce it employing the above scalar product (|l6|). Then it is necessary to introduce vectors A = A^,(0) and 
B = A™ +fc (mr). Using Fig.l and simple geometrical notions we can demonstrate the following relations in terms of 
these symbols 

1) < A • B >= |A| • |B| cosi?, costf = a(t), 

2) B = B|, +B ± ; 

3) B|, = |B|costf A = ^|B| a (i),|Bn| 2 = |A| 2 {a(t)} 2 , 

4) |BJ = |B|sintf=|B|{l-Hi)]} 1/2 , (22) 

where symbol | A| denotes the vector A length. Geometrical distance R(A, B) between two vectors A and B can also 
be found 

fc-i 

i?(A,B) = {|A - Bl 2 } 1 ^ = { £( Ai - Bj) 2 } 1 / 2 . 

Using the latter and taking into account Eqs. (||), (|l^) we can find 

R(A° k (0), AZ +k (t)) = {\K n+ kJt)\ 2 Y' 2 = v/2|A- +fc (*)|{l - a(t)y/ 2 . 

The equation above immediately shows that the distance is determined by the dynamics of evolut ion of correlation 
process. Owing to the property (0) the following relation limt^oo i?(A° (0)), A™ +fc (i)) = V2ka 2 , where a 2 is the 
variance can be developed. With regard to Eqn. (|2^) the correlation decay in limit t — > oo may result in complete 
annihilation of parallel component of state A™ +fc (t) vector . Then the state of the system at the moment t — > oo is 



to the annihilation of correlation and determined by the component B^. It results in the fact that in the limit of 
great t — ► oo the following relation: 

lrai AZ +u {t) = 0, lim A^ +fciX (t) = A™ +fc (t). (23) 

is immediately fulfilled in correspondence with to Bogolubov's |36| ] principle of correlation attenuation. 

From the physical point of view this fact means that TCF a(t) represents two interrelated states determined by 
creation and annihilation of correlation. Hence it follows that such consideration must be given to both processes in 
an explicit form for stochastic dynamics of random processes' correlation. 



IV. SPLITTING OF EQUATION OF VECTORS MOTION AND LIOUVILLIAN'S MATRIX 

PRESENTATION 

It is obvious from Eq.(22) that TCF a(t) is originated by projection of vector A™ +fc (i) ([l8|), where time t = mr on 
the initial vector of state A°(0) (see, for example, formula ([jl])). The following construction of projection operator 

UA™ +k (t) = A°(0) < A }f^^ ]> = A ?(0)«W ( 24 ) 



results from here . It is turn projection operator II from Eqn. (E4J) has the following properties 



| A°(0) >< A°(0 )[ n2 
< |A° fe (0)| 2 > 



n = 1 tl ' n 2 = n, p = 1 - n, p 2 = p, up = o, pit = o. (25) 



A pair of projection operators II and P are idempotent and mutually supplementary. Figure 1 shows that projector II 
projects on the direction A°(0), whereas the orthogonal operator P transfers all vectors to the orthogonal direction. 
Let us consider quasidynamic finite-difference Liouville's Eq. (Oh for the vector of fluctuations 



A A - +fc (*)=i%r)A£ +fc (t). (26) 

The vectors A™ +fe (£) generate the vector finite-dimensional space A(k) with scalar product in which (according to 
Eqn.(|2^)-(|2^)) the orthogonal projection operation J37J, is expressed by 

A(k)=A'(k) + A"(k), A™ +k (t) G A(k), 
A'(k)=UA(k), A"(k)=PA(k) = (l-U)A(k). (27) 

Operators II and P split Euclidean space A(k) into two mutually-orthogonal subspaces. This permits to split dynam- 
ical equation ( p6| ) into two equations within two mutually-supplementary subspaces [Q- Q as follows 

A A'(t\ 

' --iL 11 A'(t) + iL 12 A"(t), (28) 



At 
AA"{t) 

At 



= iL 21 A'{t)+iL 22 A"{t). (29) 



In the Eqs. above we crossout for short space elements indices A, A' and A" and matrix elements arguments Lij, 
Lij = IFLITj, II! = II, n 2 = P = 1 — II, i = 1, 2. In line with Refs. |39], [0 we write down Liouville's operator in 
matrix form 

? — I Lll Pl2 

\ L 2 i L 22 
Pn = IlLU, L 12 = ULP, 



To simplify Liouville's Eq. (28), (^9|) we exclude the irrelevant part A"(t) and construct closed Equation for relevant 
part A'(t). For this purpose let us solve Eqn.(p9|) step by step 



A {A ™ +fc (i)}" = zL 21 {A™ +fc (i)}' + iL 22 {A2 +k (*)}"• (31) 



Considering Eqn.(Q) we arrive at finite-difference solution of this equation in the following form 

AA"(t) 

K -L = T-^A'^t + T ) - A"{t)] = iL 2l A'{t) + iL 22 A"(t), (32) 

T 

A"(t + t) = {1 + i T L 22 }A"{t) + irL 21 A'(t). (33) 

Applying Eqs.@, @ we find 

A"(t + 2r) = (1 + iTL 22 ) 2 A"(t) + (1 + iTL 22 ){iTL 21 A'(t)} + {irL 2l }A'{t + r), (34) 

for m = 2 and 

A"(f + 3r) = (1 + irL 22 ) 3 A"(t) + (1 + iriaaf^rLaiA'^)} + (1 + irL 22 ){iTL 2 i A(i + r)} + {irL 2 iA(< + 2t}, 

(35) 

for m — 3, respectively. In general case we find 

m — 1 

A"{t + mr) = {1 + irL 22 } m A'(t) + ^ {1 + irL 22 } J {zL 21 A' (t + {m - 1 - j)r)} (36) 

i=o 

for the arbitrary number of m-steps. Then after the substitution of right side of Eq. (|36|) for Eq.(j2^) we obtain the 
closed finite-difference kinetic equation for the relevant parts of vectors 

A m+1 

—A'(t + mr) = iL n A'(t + mr) + iL 12 {\ + iTL 22 } m A"{t) - L 12 ^ {1 + irL 22 y tL 21 A' {t + (m - 1 - j)t). (37) 

j=o 

To simplify this equation, let us consider the idempotentity property, and then determine (0 < k < m — 1) 

A"(i) = 0, {I + iTL 22 } k A" (t) = 0. (38) 

Transfering from vectors A™ +fc in Eq.(^7|) to a scalar value of TCF a{t) by means of suitable projection we come to 
the closed finite-difference discrete equation for the initial TCF 

^ = l J° ) a(t) - tQ 2 J2 M ^r)a(t jr). (39) 

Here S7o is the general relaxation frequency whereas frequency uj^ describes the eigenspectrum of the Liouville's 
quasioperator L 



(0) _ < A°(0)LA°(0) > < A°L 12 £ 21 A°(0 ) > 

<|A°(0)P> ' <|A° fc (0)| 2 > 



Function M\(Jt) in the right side of Eq.(39) is the first order memory function 

_ < A ^»^'';f^; A °'°'> , M l( 0) . L (41) 
< A°(0)Li 2 L 2 iA°(0) > 



V. INTRODUCTION OF THE SET OF ORTHOGONAL RANDOM VARIABLES AND CONSTRUCTION 
OF INFINITE CHAIN OF FINITE-DIFFERENCE NON-MARKOV KINETIC EQUATIONS FOR 

DISCRETE MEMORY FUNCTIONS 



The discrete memory function Mi(jr) ( [4l| ) in Eq. (|29|) is in its turn the normalized TCF, evolution of which is 
defined by the deformed (compressed) Liouvillian's (L^ = L) 

L (1) = L% = L 22 = (1 - n)L(l - II) (42) 

for a new dynamical variable = iZ, 2 iA^(0). Thus, we can completely repeat for Mi(jr) the whole procedure 
within Eqs. (|2~4|)-(|4l|), and obtain the following non-Markov kinetic equation for the normalized TCF. The infinite 
chain of equations for the initial TCF and memory functions of increasing order results from multiple repetition of 
similar procedure. 

However this chain of equations can be obtained differently, i.e. much shorter and less costly. For this purpose let 
us employ the method developed earlier for the physical Hamilton systems with the continuous time in Refs. flip] , 
jyj. Moreover the lack of Hamiltonian and the time discretness must be taken into account. 

Let us remember that natural equation of motion ( p"4| ) is the finite- difference Liouville's equation 

— x(t) = iLx(t) (43) 

where Liouville's quasioperator is 

L = L(t, t) = {iT)- l {U{t + r,t)- 1}. (44) 

Succesively applying the quasioperator L to the dynamic variables A™ +fc (i) (t = mr, where r is a discrete time step) 
we obtain the infinite set of dynamic functions 

B„(0) = {L} n A°(0), n>l. (45) 

Using variables B„(0) one can find the formal solution of evolution Eq.(43) in the form 

A™ +fe (mr) = {1 + zrL} m A°(0) = £ ^ > .., B^O). (46) 

However, a similar form of dynamic variables is deficient. That is why we prefer the use the orthogonal variables as 
vectors W ra given below. Employing Gram-Schmidt orthogonalization procedure [42] for the set of variables B„(0) 
one can obtain the new infinite set of dynamical orthogonal variables, i.e. vectors W„ 

< W*(0),W m (0) >= 6 n , m < |W„(0)| 2 >, (47) 

where the mean < ■ • • > should be read in terms of Eqs. (|l6|)-(|l8|) and S n m is Kronecker's symbol. Now we may 
easily introduce the recurrence formula in which the senior values W„ = W„(t) are connected with the juniour values 

W o = A£(0), W a ={L-4 0) }W , W n = {L-4 n_1) }W„- 1 -^_ 1 ^ n _ 2j n > 1. (48) 

Here we used the equation, given earlier in (40) for number n=0 

Jn) < W»LW„ > <|W ?1 | 2 > 

"° - <|W W P> " n »-<|W n _ 1 |2>' (49) 

where f2 ra is the general relaxation frequency, and frequency completely describes the eigen spectrum of Liouville's 
quasioperator L (Q) . Now the arbitrary variables W n may be expressed directly through the inital variable Wo = 
Ag(0) by means of Eqs. (48) 

I L - J® ] fii ... 



The physical sense of W„ variables (vectors of state) can be cleared up in the following way. For example, in 
the continuous matter physics, the local density fluctuations may be considered as initial variables. So the local flow 
density, energy density and energy flow density fluctuations are the dynamic variables W„ where numbers n > 1. The 
careful usage of the abovementioned variables within the long-wave limits creates the basis for the condensed matter 
theory in hydrodynamic approximation. The set of the orthogonal variables (48) [see also (47)] can be connected 
with the set of projection operators. The later projects the arbitrary dynamic variable (i.e., vector of state) Y on the 
corresponding vector of the set 

n - l w " >< w »l n 2 -n p = 1 - n p 2 - p np-o 

n^Ilm ^n.m-ff-n: -^n-^m 3 nm P n , P n Ff n 0. (51) 

Let us take into consideration the fact that both sets (45) and (50) are infinite. If we execute the operations in the 
Euclidean space of dynamic variables then the formal expressions (51) must be understood as follows 

< W* Y > < YW„ > 

TT Y — W " YTT - W* " fRO") 

ll„Y-W n<|Wn|2>) Yll„- W n<|Wn|2> . (52) 

Now according to Eqs. (28)-(30), (51), (52) we can introduce the following notation for the splitting of the Liouville's 
quasioperator into the diagonal (p|™') and non-diagonal (L^) matrix elements with i =/= j, n > 1 



L(")=P„_ 1 P(- 1 )P„_ 1 , L^L^^^Luf-^, i,j = 1,2, 



n^ = n n , n 2 Il) = p n = 1 - n„. (53) 

For example, we come to the following Eqs. 

L22 = Lo = L, L22' = Pn-iPn-2 • • • P0LP0 . . . P n -iP n -\- (54) 

for the second diagonal matrix elements. Successively applying projection operators TI„ and P„ for the discrete 
equation (43) in the set of normalized TCF (t = mr) 

Mn[t) ~ < |W„(0)| 2 > (55) 

we obtain the infinite hierarchy of connected non-Markov finite- difference kinetic equations (t — mr) 

m—l 



= tJ n) M n (t) - rtf n+1 jr M n+1 (jr)M n (t jr), (56) 

3=0 

where lu^ is the eigen and fl n is the general relaxation frequency as follows 



(n) _<W*L n W n > _ („) 2 _ <\W n \ 2 > 



" <\W n \*> ' "~ 22 ' "~ <|W„-i| 2 >" 
A set of functions M n (t) (55), (56) except n = 

M (t)=a(t) = [A o (0) |2 > t = mT 

can be considered as functions characterizing the statistical memory of time correlation in the complex systems with 
discrete time. The initial TCF a(t) and the set of discrete memory functions M n (t) in Eq. (56) are of crucial role for 
the further consideration. It is convenient to rewrite the set of discrete kinetic Eqs. (56) as the infinite chain of coupled 
non-Markov discrete equations of nonlinear type for the initial discrete TCF a(t) ( discrete time t = mr everywhere) 



^ = -rn\ M x {j T )a(t - jr) + i4 0) a(t), 

3=0 

A A/f (i\ ra— 1 

^ - -rf> 2 M 2 (jr)M 1 {t - jr) + iJ^M^t), 



At 



These finite-difference Eqs. ( |5^ ) and (571) are very similar to famous Zwanzig'-Mori's chain (ZMC) of kinetic equations 
p3| - p5| , which plays the fundamental role in modern statistical physics of nonequilibrium phenomena with the smooth 
time. It should be noted that ZMC'c is true only for the physical quantum and classical systems with smooth time 
governed by Hamiltonian. Our finite-difference kinetic equations (|5^), (57) are valid for complex systems lacking 
Hamiltonian, the time being discrete and the exact equations of motion being absent. However, the "dynamics" 
and "motion" in the real complex systems are undoubtedly abundant and are immediately registered during the 
experiment. 

The first three of those Eqs.(|57j) in the whole infinite chain (|5^) form the basis for the quasihydrodynamic description 
of random processes in complex systems. 



VI. A PSEUDOHYDRODYNAMIC DESCRIPTION OF RANDOM PROCESSES IN COMPLEX 

SYSTEMS 

At first let's find the matrix elements Ly of complex systems Liouvillian's quasioperator. Employing Eqs. (0), 
(p5|),(|30|),(p3|) and ( ]54| ) we successively found 

^(0) = n a(r)- O (0) = 

T 

iL ( 2i = {r^[U(t + t, t) - 1] - a'(0)}II, (59) 
iL<f 2 ] =Tl{T- 1 [U(t + T,t)-l]-a'{0)}, (60) 
i£<°> =it- + + L®} = r- l [U{t + r, t) - 1] - 

T- X n{u{t + r,t) - 1} - r~i{u(t + r,t)- i}n + a'(o)n. (6i) 

A diagonal matrix element Lf§2 is the part of "compressed" evolution quasioperator, which in its turn is equal to 

1 + irL 22 = U{t + T,t) + ro'(0)n - {II, U(t + t, i) - 1}+, (62) 

where the anticommutator of appropriate operator is designate by the brackets {A, B} + = AB + £?AOne can see 
from the Eq.(62) that the "compressed" evolution operator differs from the natural operator U(t + T,t) because of 
the presence of contributions, associated with the first and the following derivatives of TCF the initial TCF a(t). 

The large-scale presentation of the memory function M\{t) is suitable mostly for practical applications. Using 
Eqs. (|5^)-(|6l|), and (ffi"|), ( |5~i| ) we also find the succession of the first five points of discrete functions Mi(jr) where 
j = 1,2,3,4 and 

Mi(0) = l,Mi(r) = {l/o"(0)}{o"(r) - 2ra'(0)a"(0)}. (63) 
The " Gaussian" behavior of TCF at the zero point t = 

_ <Af (0){t/ T -l}Ai°)(0)> _ 
a[U) - <|A°(0)P> 

should be taken into account in the subsequent discussion. It is proved accurately and connected directly with the 
orthogonality property of dynamical variables (^7|)-(p0|). It gives us an opportinity to simplify the memory function 
formula as follows 

Mx(0) = l, Ml (r)={^}, 
Mr(2r) = {J^}{a"(2r) - 2ra"(0)o'(r) + r[a"(0)] 2 }, 

A/ i( 3r ) = {^77TT}{a"(3T) - r a '(2r)a"(0) - 2ra"(r) a '(r) + r 2 a"(r)a"(0) + ra'(r)a"(0)}, 
a" 



Further presentation the following values M^jr) (numbers j > 5) constitute the extremely complicated combinatory 
problem. As analysis of Eqs.(^J) shows second derivative's behaviour 

Mi(jr) {-J—} a "(jr) (65) 

contributes mainly into functions M\{jr). 

Now let us move to practical realization of Eqs.(|57|), for ming a basis of pseudohydrodynamic description of corre- 
lation dynamics. Thus using orthogonal dynamic variables (|47j), ( 48 ) , (|50|) , we immediately obtain 



W = A° k , Wi={L- 4 0) }Wo = LW = (iry^Ur - l)A° k (0), 

w 2 = LWi - nlw = {l 2 - nj}Wo = (iTy 2 {u T - i} 2 A k - nlA° k , 

W 3 = LW 2 - njWi = L(L 2 - Qf)W ~ Vl\LWq = 
{L 3 - (Q 2 + n 2 )L}W = {{irf[U T - l] 3 - (ir)" 1 ^ + flj)(U T - l)}A° k . (66) 

Simple relation for the eigen and general relaxation frequencies 

4 n) - < > = o, nl= <m2> .n; = |aW(o)|, 
<\w n \ 2 > <|^„_ 1 | 2 >' 



2 _«W[0)-W o2 a( 6 )(0) - 2a( 4 )(0)(» 2 + - (» 2 + Q 2 ) 2 a^(0) 

"2- 1.9/nM > "3- a(4)(0)-(a(2)(0)) 2 ^ 



should be taken into consideration here. The orthogonal variables W n in Eq.(|6q) can be easily rearranged as follows 



Wo = A° k , W 1 =-i-A° k = -i—A k , 



= -i* 
'dt k At 



"'■< ••{■^-<«;~l!j>ijAJ + S!1)|4a°. (68) 



Those formulas ( |68| ) have considerable utility inasmuch as they permit to see the structure of formation of orthogonal 
variables and junior orders memory functions for the numbers n — 1, 2, 3. Equations (|64[), (|67|), ( |68| ) open up new fields 
of construction of quasikinetic description of random processes {AjJ(O), A™ +fc (mr)}. By analogy with hydrodynamics 

the variables Wo, W\, W 2 and W3 in Eq. ( |68| ) play the role similar to that of the local density, local flow, local energy 
density and energy flow. It is clear that this is only formal analogy and the variables W n dont possess any physical 
sense. However, such analogies can be helpful in revealing of the real sense of orthogonal variables. 

To describe pseudohydrodynamics we have to use the set of first three discrete kinetic Eqs. ( p7|) with frequences 
O? (i — 1, 2, 3) derived from Eqs. (|67|) . It is essential that all frequences fi 2 are connected straightly with the properties 
of the initial TCF a(t) only. The latter can be easily derived directly from the experimental data J45|]- |47|| . Thus 
the system of Eqs. ( pl\ ) has considerable utility for the experimental investigations of statistical memory effects and 
non-Markov processes in complex systems. 

Among them it seems to us that one could propose more physical interpretation of the different terms in the 
right side of the three Eqs. (68). For example, term — iAA/At is like a dissipation, A 2 A/ At 2 is like a inertia and 
£l 2 A(t) is like a restoring force. Third derivative A 3 Aj At 3 is the finite-difference generic form of the Abraham-Lorenz 
force corresponding to dissipation feedback due to radiative losses [see for instance formula (3) in |4g] for a recent 
experimental evidence in frictional systems]. 



VII. SHANNON ENTROPY FOR THE TIME CORRELATION AND TIME MEMORY IN COMPLEX 

SYSTEMS 



correlation . In a line with Shannon pj^ in case of discrete source of information we were able to determine a definite 
rate of generating information, namely the entropy of the underlying stochastic information by introduction fidelity 
evaluation function v{P(x 1 y)). Here the function P(x,y) is the two-dimensional distribution of random variables 
(x, y) and 

v{ p {x,y))=J J dxdyP(x,y)p{x,y), (69) 

where the function p(x, y) has the general nature of the "distance" between x and y. As pointed by Shannon [17] the 
function p(x, y) is not a "metric" in the strict sense, however, since in general it does not satisfy either p(x , y) = p(y, x) 



or p(x,y) + p{y,z) > p(x,z).. It measures how undesirable it is according to our fidelity criterion (69) to receive y 
when x transmitted. According to Shannon |l7j| any evalution of fidelity must correspond mathematically to the 
operation of a simple ordering of systems by the transmission of a signals within the certain tolerance. According to 
Shannon [fl7f the following is simple example of fidelity evaluation function 

v(P(x,y))=<(x(t)-y(t)) 2 > . (70) 

In our case it is convenient to consider the initial vector A°(0) as a variable x and the final vector A™ +k (t) at time 
t = mr for a variable y. The distance function p(x,y) [17] 

p(x : y) = ^J^ dt{x{t) - y{t)\ 2 (71) 

is the most commonly used measure of fidelity. 

Taking into accout Eqs. (|69|) , ( ffll) and the results in Section 3 as the fidelity function one can use the following 
function of geometrical distance 

v(P(A° k (0),A™ +k (t)) = 2ka 2 {l-a(t)}, (72) 

where distance function is 

p(A°(0),A™ +fc (*)) = i? 2 (A° fe (0), A" +fe (<)). (73) 



According to 1 7 partial solution of the general maximizing problem for determining the rate of generating information 



of a source can be given using Lagrange's method and considering the following functional 

P(x, y) 

{P(x, y)log pTXprjj + ^ p {x, y)p{x, y) + v{x)P{x, y)}dxdy, (74) 

where the function v{x) and \i are unknown. The following equation for the conditional probability can be obtained 
by variation on P(x, y) 



P(x,y) 
P{y) 



P v(x) = 4?f7r = ex P {-Ap(x, „)}. (75) 



This shows that with best encoding the conditional probability of a certain cause for various received y, P y (x) will 
decline exponentially with the distance function p(x, y) between values the x and y in problem. Unknown constant 
A is defined by the required fidelity, and function B{x) in the case of continuous variables obeys the normalization 
condition 

B{x) exp{— Xp(x, y)}dx = 1. (76) 

Since the distance function p(x, y) ( |7l| ) is dependent only on the vectors difference p(x, y) = p(x — y), we get a simple 
solution for the special case B(x) = a 



P y (x) = aexp{— Xp(x — y)} — aexp{— c[l — a(i)]} 



(77) 



Now let us employ Shannon's solution for continuous variables (|7q), (|77j ) a nd pass to simplified discrete two-level 
description of the system. Then let us consider the conditional probability (|77|) which describes the state on time axis 
at the moment t = hit as corresponding to the creation of correlation. Whereas the other state at the fixed moment 
t = rm which accounts for the state with the absence (annihilation) of correlation will exist. Let us introduce two 
probabilities (see Fig. 2), which will fit normalizing condition 



P 1 (i)+P 2 (i)=l, P X {t)=P cc {t), P 2 {t) = P ac {t), Pcc(t)+Pac(t) = l. (78) 

In the case of two lewels Shannon entropy 



S= -J^Pi In Pi (79) 



i=l 

increases at full disorder and takes its limiting value 

lim S = lim S(t) = In 2. (80) 

t — >oo t — >oo 

To find unknown parameters a and c in two-level description (creation and annihilation of correlation) in Eq. (|75| ) 
we should take into account normalization condition, principle of entropy increase ( |SC)| ) at t — > oo and of entropy 
extremality (presence of minimum) at full order when the following relationship: limt^ a(t) = 1 is true for the TCF. 
We obtained the following equation 

limS(t) = -{cdna + (l-a)ln(l-a)} = 
for the parameters a and c (c > 0, < a < 1) having regard to these requirements. Among two solutions (a± — 



1, ui = 0) only the first one [a.\ = 1) has physical sense. Two probabilities calculated by means of Eq.(77) will satisfy 
conditions ©, © 

Pi(t)=P cc (t) =exp{-ln2[l-a(t)]}, (81) 



P2(t) = P ac (t) =l-exp{-ln2[l-a(t)]}. (82) 

respectively. In acordance with two -lewel description it would be convenient to deal with two dynamic channels of 
entropy (creation (cc) and annihilation (ac)) of correlation (see Fig. 2) 

S cc (t) = ln2{l-a(i)}exp{-ln2[l-a(i)]}, (83) 



S ac {t) = -{l-exp[(-ln2(l-a(i))]}ln{l-exp[(-m2(l-a(t))]}. (84) 

The probabilities obtained are in the line with full dynamic (time dependent) information Shannon entropy 

S (t) = S cc (t) + S ac (t) = ln2{l - a(t)}exp{-ln2[l - a{t)}} 
-{1 - exp[(- ln2(l - a(t))]}ln{l - exp[(- In 2(1 - a(<))]}. (85) 

The entropy introduced in to Eqs. (jS3|)-(|S5|) characterized a quantitative measure of disorder in the system related 
to creation and annihilation of dynamic correlation. The probabilistic and entropy channels (|8l|)-(|84]) possess the 
following asymptotic behavior 
if a(t) -> 1, 

Pi(t) = Pcc(t) = 1 + In 2 [a(i) - l],P 2 (t) = P ac (t) = ln2[l - a{t)], 

S cc (t) S -{1 + ln2[a(t) - l]}ln{l+ln2[a(t) - l]},S ae (t) = -ln2[l - a{t)\ ln{ln2[l - a(t)}}, 
So(t) = -{1 + hi2[a(<) - l]}ln{l + ]n2[a(t) - 1]} - ln2[l - a(i)] ln{ln2[l - a(t)]}; 



Pi(t) = P cc (t) = i{l + ln2 ■ a[t)}, P 2 {t) = P ac {t) = 1{1 - In 2 • a(t)}, 

S cc (t) = ~{1 + In 2 • a(t)}ln{i[l + In 2 • a(t)}}, S ac (t) = ~{1 -ln2- a(t)}ln{^[l - In 2 ■ a(t)}}, 
So(t) = ~{1 + In 2 • a(t)}ln{i[l + In 2 • a(t)}} - ±{1 - In 2 • a(i)}ln{l[l - In 2 • a(t)]}. 



(86) 



Subsequent boundary conditions result from the equations above in terms of Bogolubov's principle of attenuation of 
correlation (^) as follows: 

limP cc =l, UmP oc (t) = 0, (87) 
lim Scc(t) = 0, lim S ac {t) = 0, lim 5 (t) = 0, (88) 

t-»0 t->0 t->0 

lim P cc (t) = 1 lim P ac = i lim S cc (t) = lim S ac (t) = lim S (t) = In 2. (89) 

It is conditions ( j37| ) that give us an opportunity to present two different states associated with the creation (cc)(in 
the time moment t — 0, P cc = 1) and annihilation (ac) (at the point t — 0, P ac (0) = 0) of correlation. Owing to 
discretness of the TCF a(t) all functions P a p, S a p as well as So(t) (a = a, c; /? = c) are discrete in the real complex 
systems. 

The results obtained in Section VI permit us to present the set of entropies for the states connected with the 
set of orthogonal variables W, and set of memory functions Afj(i) = {Mi(t), M 2 (t), M 3 (t), ■ ■ •}. On the analogy to 
Eqs.([3l|)-([S9|) these functions describe non-Markov and memory effects in the system under discussion 

Pf' (t) = P cMi (t) = cxp{ - In 2 [1 - Mi (*)]}, (90) 

P 2 Ml (t) = P Qm , (t) = 1 - exp{~ ln2[l - Afi(t)]}, (91) 

Si{t) = ln2[l - Mi(t)\ cxp{ln2[l - Mi(t)\} - 
{1 - cxp[-ln2(l - Mi)}} ln{l - exp[-ln2(l - Mi)]}, (92) 

where i — 1,2,3. Four corresponding entropies So(t), Si(t), S^i) and Ss(t) and their power frequency spectra are 
available from the set of four time functions (TCF a(t) and three memory functions Mi(t), M2(t), M^(t)). Equations 
([Sl|)-(|92"|) are of great value because they allow us to estimate stochastic dynamics of the real complex systems with 
discrete time. As a matter of principle the first three memory functions Mi(t) (i = 1,2,3) are easy to find via Eqs. 
(p7|). Using dimensionless parameter e\ = r 2 Vt\ and solution of the first finite-difference Eq.([57])we can calculate the 
discrete function Afi(jr) at the points j = 0, 1, 2, • • ■ as follows: 

Mi(0) = 1, Mi(r) = -a(2r) + e^ 1 {a{2T) - a(3r)}, 
M 1 {2t) = -{a(2T)M 1 ( T ) + a(3r)} + e^ 1 {a(3r) - a(4r)}, 
M!(3r) = -{a(2r)A//i(2r) + a(3r)Aii(r) + a(4r)Afi(0)} + e^W 4 -!-) - a(5r)}, 



Afi (tot) = - J] Mi(jr)a{(m + 1 - j)r} + e^ 1 [a{{m + 1)t} - a{(m + 2)t}]. (93) 



In general case solving the chain of Eqs. (p5|),(|57|) we can found the recurrence relations between the memory functions 
of junior and higher orders in the following form 



m — 1 



The relations obtained allow us to derive straightly the necessary memory functions M s (t) of any order s = 1, 2, ... 
from experimental data using the registered TCF a (tot) jf6|, |47). Relaxation frequencies flf,i — 1, 2, 3, given in 
Eqs. (|94| ) are available to experimental registration. Thus, it is fair to say that the applications of Eqs.(p4|) will open 
up fresh opportunities for detailed study of statistical properties of correlations in the complex systems. The very 
fact of existence of finite -difference Eqs. (||j), (|57]) enables us to develop any functions directly from the experiment. 
Therefore, the availability of discretness permits to enhance substantially the capability to get information for the 
complex systems' state. 

In conclusion let us show the equations, which characterize the rate of entropy production. It is obvious from 
conditions (|87])-(89) as well as Eqs. ([Sll)-(p5|) that the rate of entropy growth 4j| within the interval (0, oo) takes 



different sign values and is determined by the entropy behavior in the channels of creation and annihilation of 
correlation 



dS _ (ds[ 0) \ (dsi 0) \ _ fdS cc (t)\ (dS ac {i) 



dt \ dt \ dt V dt 7 V dt 



(95) 



dt 



d^^lt) 

1 K> = -ln2 V(t)exp{-ln2[l-a(t)]}{l-ln2[l-a(i)]}, (96) 



dt 

dS (t) 



\n2a'(t) exp{- La2[l - a(i)]}{l + ln[l - exp[- In 2(1 - a(t)]}, (97) 



dt 



= ln2a'(f) exp{- ln2[l - o(t)]}{ln{l - exp[- In 2(1 - a(t))}} + ln2[l - a(t)}. (98) 



The derivatives a'(t) and S' (t) here should be read in terms of Eqs. (Q), (|T3|). Since the derivative a'(t) is finite within 
the whole time interval (0, oo): |ot'(t)| < c, (where c is positive constant) the rate of entropy growth obeys the following 
boundary conditions 

^( 9 4)=0, Hmf^=0. (99) 



t-o V dt J ' t^oo y dt , 

Formulas (|95|)-(|99|) are useful for the discussion of the experimental data. Close inspection of these equations shows 
that the behaviour of derivative (^§f-) is described in many respects by the function a'(t) = T~ 1 \a(t + T) — a(t)], which 
is in its turn can be obtained from the time series observed. Relations analogous to Eqs. (j9^)- (J9^) are easily available 
for the sequence of memory functions Mj(t) (55) as well. 



VIII. APPLICATION ON ANALYSIS OF STOCHASTIC DYNAMICS OF R-R INTERVALS IN HUMAN 

ECG'S 

Let us use the stochastic dynamics of RR-intcrvals from human ECG's to illustrate of some practical value of the 
approach developed. It is well known pifl , |po| , [j49"f - |60t| that the statistical analysis of related dynamics allows the 
reliable quantitative characteristics of the human cardiovascular system states and trusty diagnostics of the various 
heart diseases |]6l| - Q. 

Most investigators into heart rate dynamics have emphasized continuous functions, whereas the heart beat itself 
is in a crucial respect a discrete event. We present here experimental evidence that by considering this quality, the 
behavior of RR intervals may be appreciated as a result of discrete dynamics. To demonstrate effectiveness of non- 
Markovian approach we only take four typical particular cases from the whole the set of experimental data |36| , which 
are available at our disposal. They are related to the case of healthy man (a), patient with a rhythm driver migration 
(b) , patient after myocardial infarction (MI) (c) , and patient after myocardial infarction (MI) with subsequent sudden 
cardial death (SCD) (d). Following standard medical practice, each from 112 person had an age, sex, and disease 
status matched pair serving as the control. 



Figs. 4-6 are figured in terms of units of r 2 . Values ei(u>) and 62(1*1) in Figs. 7, 8 are dimensionless values. Fig. 3 
shows phase trajectories, obtained for four different groups of patients in orthogonal variable's (Wo, W2) plane. Let us 
remind ourselves, that in correspondence with formulas ( |64|) , (|6^ ) the variable Wq presents RR- intervals fluctuations, 
and W2 is the second orthogonal variable and due to Eq. (p8J) is combination of an inertia force minus a restoring force. 
These variables have dimensions r and r _1 respectively, where r =< Irr > is the average value of RR-interval in time 
sequence. Set of characteristic parameters is collected in Table I. Let us mention the strong difference of numerical 
value of the first general relaxation frequency Qi frequency for four different groups of patients. Figures. 4-6 show 
power frequency spectra for three different time functions for typical patients from four different groups. Fig. 3(a) 
corresponds to a strange attractor, fig. 3(b) corresponds to quasi-periodic motion, fig. 3(c) and fig. 3(d) demonstrate 
the obviously expressed correlation of phase variables Wo and W2. Although the frequency oj is measured in units 
of and power in r 2 , respectively. Figure 4 shows power spectrum of TCF fluctuations of RR- intervals. The data, 
shown on Figs. 5, 6 are correspondingly related to power spectra of first and second memory functions. Functions 



themselves are calculated from formulas (57), ( |68| ) and (|94|). 

Figures. 7, 8 require special explanation. They show frequency spectra of first two points 6±(u>) and 62(01) of 
statistical spectra of non-Markoyity parameter (NMP) e;, where i = 1,2, ... . A presentation of NMP spectrum was 
introduced earlier in Refs. |67j , J(|| and was then used in statistical physics of liquids |^] , j70| . Close to that given 
in |37|, |6^] definitions of non-Markovity were developed later in Refs. |7l]]- [Q. In comparision with Refs. |^]- |6^] 
here we generalize NMP conception for frequency dependent case 

1 



6%(u>) = 



m(uj) 



where i = 1, 2, .. and ^i(oi) is power frequency spectrum of i-th level. 

As is shown by Yulmetyev et al in articles |^]- |7^] NMP value of 6i allows to obtain quantitative estimate of 
non-Markovity effects and statistical collective memory in random changes of experimentally measured data. 

Parameter 6i allows to divide all processes in three important cases |67j- |70|| . Markovian processes correspond 
to e > 1, while quasimarkovian processes correspond to situation with e > 1. The limit case e ~ 1 describes non- 
Markovian processes. In this case the time scale of memory processes and correlations (or junior and senior memory 
functions) coincide with each other. 

From Figs. 3-8 one can easily obtain sharp differences between four groups of patients for all types of frequency 
spectra. For instance, frequency spectrum of TCF power for healthy [Fig. 4a] is almost reproduced in NMP 6\(uj) 
spectrum given in Fig. 7a. Also it is slightly deformed in the spectra of first [Fig. 5a] and second [Fig. 6a] memory 
functions and is strongly transformed in NMP 62(10) spectrum (Fig. 8a). Sharp peak in the vicinity of the point 
with u> ~ 0.125/.U, being characteristical for the patient (b), is seen in the power spectrum of first and second MF's 
(Fig. 5(b), 6(b)). However for other spectra of type b (for example Figs. 6b, 7b, 8b) quite complicated structure 
appears. Frequency spectrum of type (c), which is characteristic for IM, contains two sharply expressed spectral peaks 
nearly the frequencies, approximately 0.2 and 0.4 f.u. on the background of low intensity white noise . These peaks 
are conserved in the spectra of first (Fig. 5c) and second (Fig. 6c) MF. In NMP spectra 6i(oo), 62(00), complicated 
structure of spectral lines also appears. In characteristic case of patient with SCD frequency spectra of type (d) 
everywhere contain sharp peaks close to frequency 0.25 f.u. We would like to mention that all frequency spectra (5,6,7 
and 8) are persuasively saying for strongly expressed non-Markovity for time change of RR- intervals. 

Figs 7a-d and 8a-d shows, that all values of NMP 6\(ui) and 62(01) lie in small interval of values (0 4- 30). 

This fact convincingly tell us about characteristic statistic memory and noticeable non-Markovity effects in statistical 
dynamics of RR-intervals from human ECG's. Obtained results on non-Markovian properties of temporal behavior of 
RR-intervals justify significant and characteristic differences in data for all 4 groups of patients. We hope that the use 
of non-Markovian dynamics in the spirit of developed theory will incorporate development of more precise estimate 
of the state of cardiovascular systems for healthy as well as for more careful diagnostics of different patients. 



IX. DISCUSSION 



The present paper deals with two interrelated important results. The first one is connected with the establishment 
of the chain of finite-difference non-Markov kinetic equations for the discrete TCF. In this case the state of complex 
systems at the definite level of correlation is described by two vectors constructed over the strict determined rules. It 



operations and matrix elements of Liouvillc's quasiopcrator ensure the splitting of natural equations of motion and 
then they are solved in the closed finite-difference form. Using Gram-Schmidt orthogonalization procedure we find 
an infinite set of the orthogonal dynamic random variables. This it allows us to obtain the whole infinite chain of 
finite-difference kinetic equations for the initial discrete TCF. These equations contain the set of all memory functions 
characterizing the complete spectrum of non-Markov processes and statistical memory effects in the complex system. 
The presence of discretness and the very fact of the existence of finite-difference structure enable, in principle, to find 
all memory functions solving successively kinetic equations for the TCF. Parameters of these equations can be easily 
obtained from the experimentally registered TCF. In chaotic dynamics of complex systems the TCF above plays the 
role similar to that of the statistical integral in equilibrium statistical physics. 

Another important result of our work is the dynamic (time dependent) information Shannon entropy given in 
terms of the TCF. It allows us to use the information measure for the quantitative characteristic of two interrelated 
correlation channels. One of them corresponds to the creation of time correlation and the other to the annihilation 
of correlation. 

For that as we employ one of the classical Shannon's results 17 1, related to the introduction of fidelity evalution 
function and distance function between two vectors of state. The existence of a new information measure opens up 
new fields for exploration of information characteristics of complex systems. In particular, some interesting data arise 
from calculations frequency spectra of power of information entropy. 

The important consequence of the results obtained is the usage of power spectra of memory functions Mj(mr), 
where m — 0, 1, 2, 3, ... and j = 1, 2, 3, .... The set of three junior memory functions with numbers j = 1 , 2, 3 provides 
the basis for the pseudohydrodynamical description of the complex system. In practice, any memory function can be 
extracted from the experimental time sets and experimentally recorded TCF. These criteria provide the possibility to 
get reliable information about non-Markov processes and memory effects in natural evolution of complex systems. In 
principle, the new point in the analysis of complex systems arises from the opportunity to construct the dynamical 
information Shannon entropy for the experimental memory functions. Undoubtedly, detection of the frequency spectra 
of power of entropy for memory functions gives us new unique information about the statistical non-Markov properties 
as well as memory effects in complex systems of various nature. 

Application of the theory developed on the analysis of dynamics of RRintervals from human ECG's strongly suggest 
the substantially non-Markovian properties of the this dynamics. Here we have obtained non-Markovian quantitative 
characteristics for the fourth various groups of patients. One might expect this method may be use in distinguishing 
healthy from pathalogic data sets based in differences in its non-Markovian properties. 

In conclusion it may be said that this paper describes a first-principle derivation of a hierarchy of finite-difference 
equations for time correlation function of out-of-equilibrium systems without Hamiltonian. The approach developed 
seems to have potentials and offer few advantages over the usual Hamiltonian point of view. A similar situation are 
true apparently with regard to turbulence, aging, for istance, as in spin glasses and glasses as well as experimental 
time series for living, social and natural complex systems (physiology, cardiology, finance, psychology, and seismology, 
etc.) 

By way of illustration it is significant that the anomalous scaling of simultaneous correlation function in turbulence 
is intimately related to the breaking of temporal scale invariance, which is equivalent to the appearence of infinitely 
many time times scales in the time dependence of time-correlation functions. In Refs. ]75|| was addresed temporal 



multiscaling on the basis of the continued fraction representation of turbulent correlation function 76 within the 
framework the Zwanzig-Mori formalism fl43| , [ H w hich was applied to the time correlation function in turbulence. 
It has been shown by Grossman and Thomas |76|] that the Zwanzig-Mori formalism applied to turbulent systems 
described by Navier- Stokes-like equations. 

Mode coupling equatins have been considered in various areas of many particle physics for an approximate treatment 
of the dynamics of particles in glasses |7^] , |78| . These equations are obtained if one represent within the Zwanzig-Mori 
formalism |43| , correlation functions in terms of memory kernels and then expressed the latter via a factorization 
approximation in terms of the former for the glass transition of molecular liquids ]79| ] . It was been found by Heuer 
et al. |8(]] that a model-free interpretation of higer-order correlation function determined by NMR reveals important 
information about the complex dynamics close to glass transition of polymers. This is now been demonstrated with 
spin glasses that how a hierarchical model of spin glasses relaxation can display aging behavior in the time scale, 
similarly to what is found in spin glasses and other complex systems out of thermodynamical equilibrium. 

The application of the approach developed on the analysis of the temporal behavior of complex systems of various 
nature will be available in our forthcoming papers. 
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XI. FIGURE CAPTIONS 



Fig.l. Simple geometrical notion on vectors, their scalar product and normalized TCF of random variables. 

Fig. 2. Scheme of simplified two-lewel description of a complex systems state. Two probability P cc {t) and P a c(t) 
decsribes a stochastic processes of creation (existence) and annihilation (decay) of time correlation. 

Fig 3. Phase-time portreit in orthogonal variables (Wo, W%) plain (see formulas (|66|), ( |6S| ) for fourth group of 
patients: healthy (a), patient with rhythm driver migration (b), patient after myocardial infarction (c) , and patient 
after MI with subsequent SCD (d). As a matter of fact us utilised dimensionless variables Wq/t and Wa/r -1 . 

Fig 4. Frequency spectrum of power [Iq{uj) for TCF of fluctuation of RR- intervals for fourth patient groups: 
healthy (a), patient with rhythm driver migration (b), patient after myocardial infarction (c) , and patient after MI 
with subsequent SCD (d). The schedule is submitted in dimensionless units. The frequency is marked in terms of 
units of (2tt/t), the function fj,o(uj) is figured in units of r 2 . 

Fig 5. Frequency spectrum of power fii((J) for the first MF M\(t) for fourth patient groups: healthy (a), patient 
with rhythm driver migration (b), patient after myocardial infarction (c) , and patient after MI with subsequent SCD 
(d). The schedule is submitted in dimensionless units. The frequency is marked in terms of units of (2tt/t) , the 
function is figured in units of r 2 . 

Fig 6. Frequency spectrum of power Hi(oS) for the second MF M 2 {t) for fourth patient groups:healthy (a), patient 
with rhythm driver migration (b), patient after myocardial infarction (c) , and patient after MI with subsequent SCD 
(d). The schedule is submitted in dimensionless units. The frequency is marked in terms of units of (2tt/t) , the 
function /^(w) is figured in units of r 2 . 

Fig 7. Frequency spectrum of the first point in the statistical spectrum on non-Markovity parameter ex(u>) for 
fourth patient groups:hcalthy (a), patient with rhythm driver migration (b), patient after myocardial infarction (c) , 
and patient after MI with subsequent SCD (d). The schedule is submitted in dimensionless units. The frequency is 
marked in terms of units of r 2 . 

Fig 8. Frequency spectrum of the second point in the statistical spectrum on non-Markovity parameter £2(0;) for 
fourth patient groups:healthy (a), patient with rhythm driver migration (b), patient after myocardial infarction (c) , 
and patient after MI with subsequent SCD (d). The schedule is submitted in dimensionless units. The frequency is 
marked in terms of units of r 2 . 
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